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Abstract 

The question of the uniqueness of the Brodsky-Lepage-Mackenzie procedure for fixing the renor- 
malization scale in perturbative QCD is discussed. It is shown that the resulting finite order 
approximants are as ambiguous as the original truncated perturbative expansions. This inherent 
ambiguity of the BLM procedure undermines the recent attempts to define "genuine" higher order 
perturbative corrections. 



1 Introduction 

Over the last 15 years the problem of the renormalization prescription dependence of finite order 
approximants to perturbation expansions in QCD Q has been a subject of lively and sometimes 
even heated debate. From time to time a "resolution" of this problem is announced, but invariably 
it turns out that these "solutions" contain the original ambiguity in some disguise. One of such 
methods has been suggested some years ago by Brodsky, Lepage and Mackenzie (BLM) [IJ. Soon 
after the appearance of the preprint version of [|J Stevenson and Celmaster Q pointed out its 
inherent ambiguity, but their criticism has largely been ignored and the BLM procedure has since 
been used in many phenomenological analyses. In the journal version of their work the authors of 
[|l|] have acknowledged the presence of the ambiguity pointed out in [Q] , but claimed that it can "be 
eliminated to a large extent by adopting some physical process as a theoretical standard for defining 
a s (Q)." As the BLM procedure is based on the generalization (called "naive nonabelianization" 
]|]) of the QED procedure of incorporating the effects of quark loops in the renormalized coupling 
constant (couplant for short), it has been conjectured to be closely related to the large order 
behaviour of the expansion coefficients in perturbation theory. This observation has recently been 
exploited in several attempts ||, [|, |5) to help answer the question of the importance of higher order 
corrections in QCD perturbative expansions. The idea suggested in these papers is to distinguish 
the "genuine" higher order corrections from those governed by the choice of the renormalization 
scale and, supposedly, related to the large order behaviour of perturbative expansions. As this 
question is of principal as well as practical importance, it is certainly useful to reappraise the old 
criticism voiced in [Q] in order to see how much of it is relevant for these recent efforts. The main 
purpose of this paper is to carry out such an analysis. 

In this article only the renormalization prescription dependence of physical quantities in QCD with n / massless 
quark flavours will be considered. 



The paper is organized as follows. After introducing in Section 2 the notation and recalling some 
basic facts concerning the description of ambiguities of finite order perturbative approximants, I 
shall discuss in Section 3 the BLM procedure in QCD. This is followed in Section 4 by a short 
detour to QED, in order to see why there the BLM procedure does, indeed, lead to a unique result. 
In Section 5 the abovementioned suggestion to avoid the "residual" RS dependence of the BLM 
results by selecting some physical process as a standard for the definition of a s is shown to be no 
real cure to the problem. The implications of the inherent ambiguity of the BLM procedure for the 
basic strategy of the papers ||, [|, |j are discussed in Section 6, followed by a short summary and 
conclusions in Section 7. 



2 Ambiguities of finite order approximants in perturbative QCD 

Before coming to the essence of the BLM procedure, let me introduce notation and recall a few 
basic facts concerning the renormalization prescription ambiguity of finite order perturbative ap- 
proximants. I shall make a clear difference between the concepts of renormalization prescription, 
renormalization scheme and renormalization convention, which all appear in the literature and 
which are not always used in the same sense. I shall emphasize the mathematically complete 
quantitative description of these concepts as this will be crucial for the later discussion. Let us 
start by considering perturbation expansion for the generic physical quantity r(Q), depending, for 
simplicity, on a single external momentum Q, 

r{Q) = a(p, RS) [l + n(Q/fi, RS)a(/j, RS) + r 2 (Q/(i, RS)a 2 (/x, RS) + • • •] . (1) 

The generalization of the following considerations for the case when the leading term behaves as 
a p is trivial. In the above expression a(/x, RS) = a s /ir is the renormalized couplant (the adjective 
"renormalized" will be dropped in the following), evaluated at the renormalization scale /j in a 
given renormalization scheme (RS), to be fully specified below. While the whole sum in ([l]) is 
independent of both the scale /i and the RS, any finite order approximant 

AT-l 

rW(Q/M, RS) = £ r fe (Q/ / x,RS)a fc+1 (^,RS), (2) 

fc=0 

does depend on /i as well as the RS. The dependence of the couplant on the scale [i is governed by 
the equation 



da(/i, RS) 
dln/x 



/3(a) = -ba 2 (fi, RS) (l + ca{fi, RS) + c 2 a 2 (/i, RS) + • • •) . (3) 



In massless QCD the first two coefficients on the r.h.s. of @ are unique functions of n/, the number 
of massless quarks (b = (33 — 2nj)/6, c = (153 — 19ny)/(33 — 4n^)), while all the higher order 
ones are completely arbitrary 0. Once they are chosen and some initial condition on a is specified, 
@ can be solved. One way of specifying this boundary condition is via the scale parameter Ars> 
introduced in the following implicit equation for the solution of (||) Q [] 



b In - — = — h c In h / dx 

A RS a 1 + cfl Jo 



+ 



x 



2 B( n ){x) x 2 (l + cx) 



(4) 



2 These coefficients define the so called renormalization convention (RC), introduced in |J. 

3 The parameter A introduced in (^J) differs from A used in most phenomenological analyses by multiplicative 
n/-dependent factor (2c/b) c ^ b , which for realistic values of nf is close to unity. 
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where B( n \x) = (I + CX + C2X 2 + ■ ■ -+c n _ix" 1 ). The dependence of the couplant on the parameters 
Q for i > 2 is determined by equations similar to (|3|) || 

M^Ci) _ bx i+2 

Pi = -f3{a) / y^ dx, (5) 



dc, ~" ^ >J (/3(x))2 

which are uniquely determined by the basic /3-function in (||) and thus introduce no additional 
ambiguity. The renormalization scheme is defined by choosing both the coefficients c, and the 
parameter A^g, i.e. by the full specification of the solution of (Q), RS= {q, A^g}. Note that at the 
NLO only A^g labels the RS. Finally the renormalization prescription is defined by the specification 
of both the renormalization scale [i and renormalization scheme RS. Only the specification of the 
renormalization prescription thus leads to unique results of finite order perturbative approximants 
©• 

In connection with Apj,g the following fact is worth mentioning. The BLM procedure, to be 
discussed in the next section, is based on the isolation of the n^-dependence of the expansion 
coefficients r^. This brings up the subtle question of the n/-dependence of the chosen RS, which 
up to the NLO means the ny- dependence of Ajj,g. As worlds with different number of massless 
quarks are not related by any theoretical arguments, such as the renormalization group (RG) 
considerations, there is, however, no meaningful way of introducing this dependence. What can, 
however, be done in a well-defined manner is to discuss the n/-dependence of the ratio A RS //ARg 
of the A-parameters, corresponding to two different RS, as this ratio is determined solely by the 
RG considerations. 

After the arbitrary coefficients Cj are specified there are thus two parameters to vary: \i and 
Aftg. As, however, [i enters the solution of (||) always in the ratio with A^g, we can 

• select one A^g and thus one of the RS={q, Ars} and vary fi only, 

• fix fi by identifying it with some external momentum, such as Q, and vary A^g instead. 

Both of these options are completely equivalent and it is merely a matter of taste which one to use. 
To vary both the scale fi and the A^g is legal, but obviously redundant. Fixing the scale \i without 
also specifying the RS is, on the other hand, insufficient to specify the renormalization prescription 
and thus does not lead to a unique result for the finite order approximants (Q). 

I emphasize this point as the scale /i is often fixed by identifying it with some natural physical 
scale of the process, such as Q. Although such a natural scale can usually be found, its mere 
existence does not fix the RS and thus does not specify the couplant. In most of phenomenological 
analyses the MS RS is tacitly adopted, but there is no theoretical argument for this choice, except 
that in this RS and for the conventional choices of \x = Q, the coefficients are often, but not 
always, small. If, however, the magnitude of the expansion coefficients should be the criterion for 
the "best" choice of the scale fi and the RS, we would be naturally drawn to the effective charges 
approach 0, where all higher order coefficients r/% vanish by definition. 

While the explicit dependence of the couplant on the renormalization scale \x and the parameters 
ARg , c% , i > 2 is given by (||) and (||) , the dependence of the coefficients on them is determined 
by the relations 

drW(QMRS) N+1 drW(Q//x,RS) N+1 drW(Q//x,RS) N+l 

din// ~ U[a h dlnA RS ~ U[a h da ~ U[a h W 

which express the internal consistency of perturbation theory. Iterating these equations we easily 
find 

n(Q//i, RS) = &ln^--p(Q/AR S ), 
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r 2 (Q/p,RS) = p 2 -c 2 + r 1 +cr 1 , (7) 

r n (Q//x,RS) = p n - c„ + f{ri,Ci,pi\i < n - 1). 

where the quantities p, pi are RG invariants. The dependence of the perturbative approximants (^) 
on Q comes entirely through the invariant p(Q/A), which can be written as 

p = 61n(Q/A RS )-ri(l,RS), (8) 

while all the higher order invariants pi are just pure numbers. A nontrivial part of any perturbative 
calculation beyond the LO boils down to the evaluation of the invariants pi, the rest being essentially 
a straightforward exploitation of the RG considerations based on (3,5,7). 

In the following section the BLM procedure will be discussed in detail at the NLO, as this is 
where the dependence on the scale p comes in and where also its basic ambiguity becomes evident. 
At this order only the first two, unique, coefficients b and c in (J3j) are taken into account and the 
renormalization prescription dependence of (]!]) is therefore essentially a one-parameter ambiguity. 
Varying either p for fixed Ars, or vice versa, we find 

a(//,RS) = a(p,RS) (^1 - bin (^J a(p,RS) + ■ ■ , (9) 

a(p,RS') = a(p,RS) (l-bhi(j^J a(/x,RS). + • • (10) 

In both ways of labelling the ambiguity, we have 

a = a (1 — na H ) , r[ = r\ + k (11) 

At the NLO and for some "initial" a = a(p, RS) n can thus be used as yet another way of labelling 
this ambiguity. After selecting 

a) the scale p in a fixed RS, or 

b) the parameter Ars for fixed p, or 



c) the parameter k for a fixed initial a in (11) 



we should, however, check that the resulting couplant and expansion coefficient r\ do not depend 
on the RS (for a)), p (for b)) or initial a (for c)). At the NLO and in all three ways of labelling 
the inherent one-parameter ambiguity, both the Principle of Minimal Sensitivity (PMS) Q and 
Effective Charges (ECH) approaches do, indeed, lead to a unique result for A 2 \Q) 0. On the 
other hand, as will be shown in detail in the next section, this basic requirement is not satisfied by 
the BLM procedure. 

Although the scale p appears naturally and unavoidably in the process of renormalization, 
perturbation expansions for physical quantities can actually do without it. Combining eqs. (||) 
and @ allows us to express all the expansion coefficients as unique functions of the /3-function 
coefficients c& and the couplant a. To specify a unique result for any finite order approximant (||) 
we can thus use the set {a, q} instead of {p, Ars, Cj}. At the NLO this means that instead of the 
pair p, Ars the couplant itself can serve to label the one-parameter ambiguity! Note that neither 
the PMS nor the ECH approaches fix the scale p and the RS={Ars} separately, but merely their 
ratio, or, equivalently, just the couplant a. 

4 Despite the appearance of Ars in this expression, p is actually RS-independent as the RS-dependences of the 
two terms on the r.h.s. of mutually cancel. 

5 Starting at the NNLO, there are certain complications with the existence and/or uniqueness of the PMS and 
ECH "optimized" results, but this is irrelevant for the present discussion. 
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3 The BLM procedure in QCD 



This method borrows its basic idea from QED, where the renormalization of the electric charge 
is entirely due to the vacuum polarization by the charged fermion-antifermion pairs. In QCD the 
scale dependence of the couplant is due to other effects as well, but the authors of [|J suggest to 
fix the scale fi by absorbing the effects of quark loops entirely in the definition of the renormalized 
couplant. Let us consider the quantity ([[]) up to the NLO and define the class of the so-called 
regular RS={Ars} by the condition that the expansion coefficient r\ can be written as a linear 
function of nj 



r (2) (Q//i,RS) =a(yu,RS) (l + r 10 Rs) +n f r n (|,RS 



a(/i,RS) , 



(12) 



where no> r ii ar e ny-independent coefficients []. The BLM procedure fixes [i by the condition 
rn(/VQ,RS) = 0. This implies § 



^blm(/^,RS) 
n(MBLM/Q,RS) 



^exp [3rn(/i/Q,RS)] , 
33 

no ( M /Q,RS) + yru (/x/Q.RS) 



(13) 
(14) 



and consequently 



rg M (Q/ M ,RS) 



°(/ i BLM) RS) ^1 + 



.-,.(§,Rs) + fr u (f,B S 



o> (mblmj RS) 



(15) 



Note that the BLM-fixed scale /iBLM is a function of both the "initial" scale and RS. This by 
itself is not surprising as the same happens in the PMS and ECH approaches as well. In these 
approaches the corresponding couplant a, as well as the coefficients r^, are, however, independent 
of both [i and the RS, while this is not the case in the BLM procedure. In the explicit calculation 
[§] Celmaster and Stevenson have shown that for the T hadronic decay width the coefficient t\ 
has different values in two most frequently used RS's, namely the MS and symmetric MOM RS 
based on the 3-gluon vertex in the Landau gauge. Despite the fact that this simple example is 
sufficient to demonstrate the basic shortcoming of the BLM procedure, its message has largely been 
ignored, presumably because the dependence of BLM results on the RS has been considered as a 
kind of "residual" dependence, less important than the dependence on the scale \x. In the rest of 
this section I shall demonstrate that this "residual" RS dependence of the BLM procedure actually 
coincides with the original scale ambiguity in a fixed RS. 

To see how the BLM results depend on \i and the RS, we need to evaluate the ratio 



/"BLM(/i,RS) 

Ars 



exp 



exp[Z(/i,RS)]. 



(16) 



Going from RS to RS' (for fixed /i), or from \x to // (for fixed RS), we get, exploiting (|7|), 



3ri(/VQ,RS') = 3r 1 ^/Q,RS)+(y-n / jln. ; . 

33 \, fn' 



33 



/Ars\ 



3n(///Q,RS) = 3ri(/VQ,RS) + ( — -n f J In ^ 



(17) 
(18) 



6 In this statement fi as well as the RS (labelled by Ars) are held fixed and only the explicit dependence of n on 
rif is taken into account. 
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Within the class of regular RS's we have 



bluffs) = A + Bnj, (19) 

Z(p,BS') = Z(m,RS) + ~ (20) 

33 

J"i (a*blm(M) RS')) = ri (/j-blm(M) RS)) + k, k = —B + A, (21) 

oblm(RS') = oblm(RS) [1 - koblm(RS)] . (22) 



The last two equations coincide with the transformations (|11|). We see that the BLM results are 
unique only provided k = 33B/2 + A = 0, which is equivalent to the condition that the ratio 
Ars/Ars' is rij-independent. This condition divides all regular RS into disjoint subclasses of RS's, 
characterized by the property that for any pair of RS's from the same subclass k = 33B/2 + A = 0, 
while for RS, RS' from different subclasses k / 0. The BLM results are thus unique only within 
each of these subclasses, while different subclasses lead to different results. To show that the latter 
situation arises naturally for quite conventional RS's, let me take RS=MS, which is the standard 
choice in all BLM papers and consider for RS' the set of MOM-based RS's, with the couplant 
defined via the three-gluon vertex at the symmetric point and for general gauge fixing parameter 
etc |§ and denote it RS(ac)- While the gauge parameter a G must also be renormalized, it can, as 
argued in 0, be considered as constant, i.e. not running. In fact only in such case [] are the two 
lowest order /3-function coefficients b, c RS-invariant constants. Within this class of MOM-based 
RS's, the corresponding Amom(«g)j and thus also the couplant and expansion coefficient ri, become 
functions of a G € (— oo,+oo). Consequently the ratio Ars/Ars' is manifestly ny-dependent and 
the expression for k in (21) can be read off eq. (18) of Q: 



Kioto) = — I ~ — (1 - I)a G -(-- + -)a 2 G - —a G (23) 
v UJ 48 16 v ; V 8 8/ G 16 G y ' 

where I = 2.3439. Varying the gauge fixing parameter etc between — oo and +oo makes k(«g) to 
span the same interval. For any given r\ one can thus always find such a value of a G that when 
k{olq) is used on the r.h.s. of (pi]), the result is just this r\. Note that at the NLO the MS RS is just 
a special case of the MOM-based RS(«g(MS)) with «g(MS) given by the solution to the equation 
k{o.q) = 0! So in any RS and for any scale /i, r( 2 )(/z, RS) coincides with the result obtained via 
MOM-based BLM procedure for some gauge parameter ac^RS): 

( 2 V/Q,RS) = r( 2 )(/x B LM/Q,RS'(a G (/x,RS))) (24) 



r 



We conclude that not only do the BLM results depend on the chosen RS(ac), as claimed in H, 
but quantitatively this ambiguity is completely equivalent to the scale ambiguity in any fixed RS. 

In analogous way it can be shown that for nj-independent scale [i the BLM results are, indeed, 
yu-independent. Formally we could introduce ny-dependent scale /i as well and thus destroy the 
scale independence of the BLM results, but as for each fixed n/ the scale \x is completely arbitrary 
there is no meaning in introducing this dependence. The situation looks differently as far as the 
n/-dependence of the chosen RS's is concerned, primarily because in the MOM-based RS's this 
dependence arises quite naturally. 

It is also simple to see why, contrary to the BLM approach, the ECH one 0| leads to unique 
results. Taking into account ([?]) the condition r\ = implies 

Mech(RS) = P(g)_ (25) 
Ars b 



7 Or for the Landau gauge ac = 0. 
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The scale /i does not appear in ( pq ) at all and as the r.h.s. of (25) contains only RG invariants 
p(Q) and b, it is manifestly RS-independent. Similar considerations hold in the PMS approach. It 
is thus the very essence of the BLM procedure, i.e. the attempt to separate the coefficient r\ into 
the nj-dependent and nj- independent parts, that causes its inherent ambiguity. 



4 The case of QED 



The results of the BLM procedure in QED can be obtained from (po|)-(21) simply by taking into 
account that in QED the corresponding 6qed = — 2n//3 and contains thus no term analogous to 
33/2 = lliV c /2, coming in QCD from gluon selfinter action. Consequently the BLM results are 
unique provided k = A = 0, which is again equivalent to the condition that Ars is n/-independent. 
In QED, however, the scale as well as RS-dependence of the corresponding couplant come entirely 
from the renormalization of fermion loops and it is therefore natural to define the regular RS's as 
those satisfying this condition. In QED and within any class of regular RS's the BLM results are, 
indeed, unique. 



5 Fixing the RS by means of the "standard" physical process 

The idea, suggested in Q, is to fix the RS of the couplant a(fi, RS) with the help of some "standard" 
physical quantity, such as the familiar ratio 

cr(e+e -»■ n+/i ) 

by demanding r(s) = a(y / s, RS). Note that in this way chosen couplant is just the effective charge 
corresponding to the physical quantity r(s). 

The main problem with this strategy is that there is no theoretical reason to prefer one particular 
physical quantity to another to serve as the "standard" . By appealing to some "standard" physical 



process to fix the RS this ambiguity is transformed into the "initial condition" ambiguity flO| , 11 1. 
Moreover, if the ECH approach is used for the "standard" physical quantity, why not to use it for 
the one under study? 

Secondly, in writing r(s) = a( v / s, RS) we have already set the scale of the couplant a(/i, RS) to 
/i = y/s. Recall that, as emphasized in Section 3, the ECH approach does not actually fix the scale 
and the RS, but determines directly the couplant oecH) or equivalently, the ratio /x/Ars! Without 
specifying the scale \x in the equation r(s) = oech(^,RS) the RS is not fixed and the reference to 



the "standard" physical quantity thus of no help. However, to set fj> = y/s in (26) has no theoretical 
justification, in particular taking into account that the basic aim of the BLM procedure is just to 
find some plausible scale fixing method! 



6 Generalization of the BLM procedure and "genuine" higher or- 
der corrections 

Let me now turn to the implications of the inherent ambiguity of the BLM procedure for some of 
the essential ingredients of the recent papers ||, |5| , where higher order perturbative corrections 
are separated into two parts: 

• the so called "genuine" higher order corrections, which are "hard to anticipate" |J|, and 
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• those related to the renormalization scale dependence of the couplant, which by an "improper" 
choice of /j, may become artificially large. 

The starting point of these attempts is the claim ||] that "some prescriptions may be closer to 
general expectations." Without specifying the meaning of the term "general expectation", all the 
mentioned papers single out the BLM procedure as the best way of fixing the scale fi. One of the 
arguments in favour of this conjecture is the widely held view that because the BLM approach 
fixes the scale [i via the "naive nonabelianization" of the original QED procedure of incorporating 
the effects of quark loops in the renormalized couplant, it is related to the large order behaviour of 
expansion coefficients in (||) || . By generalizing the BLM procedure to higher orders the papers 
[]|, [ij] attempt to incorporate the effects of the leading IR renormalon in the BLM couplant, thereby 
isolating the remaining "genuine" higher order perturbative corrections. 

The generalization of the BLM scale-fixing procedure to arbitrary an order N suggested in ||, Q 
is based on the generalization of the relation (|l~2] ) 

r n = r n o + r n in f H h r nn n n } . (27) 

Writing r n as 

r n = $n + {b/2) n d n , (28) 

where d n are rif- independent, absorbs the leading nj-dependence in d n and allows the authors of 
[||, to define the generalized BLM scale at order N via the relation f\ 

N 

a(Q { ^ M ) = a(Q)M^, M N = 1 + £ d n (b/2) n a n . (29) 

n=l 

With these definitions eq. ([[]) can be written as 

N 

r(Q) = a(Qff M ) + £ 5 n a(Q) n+1 (30) 



n=l 



Except for the argument of the couplant in the second term of (^), which, however, does not 
influence the coefficient 5\ = ri(/UBLM/Q> RS), ( |30| ) reduces for N = 1 to (15). According to [||, |2| 



the coefficients 5 n represent the "genuine" higher order corrections, in contrast to those included 
in Mjv, which are incorporated, via the BLM procedure, in the leading-order couplant o,(Q^ M ). 

Moreover, in the limit — > oo, Qblm = hniAr^oo Qblm ^ s claimed to be independent of the finite 
renormalization of the fermion loop. 

However, in QCD the renormalization of the colour charge is not given by the gluon polarization 
due to fermion loops only. Therefore the fact that each of the terms on the r.h.s. of ( pop is 
independent of the finite renormalization of the basic fermion loop does not imply that it is also 
RS-invariant and thus unique. In Section 3 I have shown explicitly that within the class of MOM- 

(2) 

based RS's both the BLM couplant o(Q B l M ) and the leading coefficient 5\ of the second term in 
(|3^) are RS-dependent. The value of the "genuine" NLO correction coefficient <5i is a function of 
the the gauge-fixing parameter ac, and can take any prescribed value. The RS-dependence of 5± 
automatically implies that the separation of r(Q) into the two terms in (|3(]) is ambiguous even if 
the BLM procedure is generalized to an arbitrary order N. Recall that in perturbative expansions 



8 As the leading order coefficient ro of Q ^] is unique, it can be set to unity, as done in (^), without losing any 
generality. 

following H [| I set in the rest of this section /i = Q. 
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of RS-independent quantities J2 n =p r n« n the leading-order coefficient r p has to be unique as there 
is no way how its RS-dependence could be compensated by the RS-dependence of the expansion 
parameter a! The only way to secure the RS-invariance of the full sum in (|3^) is therefore the 
mutual cancellation of this dependence between the BLM couplant a(QeLM) an d the expansion 
describing the supposed "genuine" higher order perturbative corrections. 

Beyond the NLO the couplant a as well as the expansion coefficients become functions of the 
additional free parameters q; i > 2, specifying the RC. This freedom implies that not only 6\, but 
in fact all the coefficients 5^ are completely arbitrary, exactly in the same way as the coefficients 
r& in the general RS! In the ECH approach, for instance, they are set to zero. The large order 
behaviour of the coefficients r^ crucially depends on the choice of the RS and so do also all the 
higher order coefficients 5^. The reason is clear: the RG transformations mix the two terms on the 
r.h.s. of ( pC)| ) so that we cannot attach an unambiguous meaning to each of them separately. 

The fact that by an appropriate choice of the RS we can get rid of the divergent behaviour 
of rk does not, of course, mean that we can in this way solve the problem of the divergence of 
perturbation expansions. The standard derivation of the asymptotic behaviour of the expansion 
coefficients is carried out in the 't Hooft RS, in which all higher order /3-function coefficients 
Cj-; i > 2 are set to zero by definition. The resulting growth behaviour of the coefficients at large 
k implies via factorial behaviour of the invariants pk oc k\. This, in turn, implies that are 
f actor ially divergent also in all those RS={Apj,g, q}, where the coefficients c& define a convergent 
series. In the RS's where r\~ define convergent series, the factorial behaviour of the invariants pk 
reappears as the divergence of the perturbation expansion of the corresponding /3-function. Note, 
however, that factorial behaviour of the coefficients Ck in (||) cannot be influenced by the choice of 
the scale p and thus has nothing to do with it. The fact that we can freely shuffle part or the whole 
factorial divergence of the RG invariants p^ between the coefficients and clearly signals that it 
is impossible to relate the effects of the factorial divergence of the invariants pk to the choice of the 
scale, as suggested in || f|, ||. In other words, there is no well-defined relation between the large 
order behaviour of the coefficients and the choice of the renormalization scale p. The former 
depends not only on p but also on all the parameters specifying the RS. 

It is also fair to say that we actually do not even know how the coefficients Ck behave in the 
most popular RS, the MS. The usual expectation is that they define a divergent series, but there 
are no convincing arguments behind this conjecture. And without this knowledge the behaviour of 
the expansion coefficients in this RS is also an open question. 

7 Summary and conclusions 

The preceding sections demonstrate two closely related facts: 

a) the inherent ambiguity of the BLM scale fixing procedure, and 

b) the impossibility to define in a reasonably unambiguous (i.e. RS-independent) way the "gen- 

uine" higher order perturbative corrections to physical quantities. 

The RG transformations bind inextricably the two terms in the sum of (^0|). The "genuine" higher 
order corrections, defined in [||, |5| by incorporating the supposed large order behaviour of the 
expansion coefficients in the BLM scale, are in fact as ambiguous as the original perturbation 
expansions before the BLM scale-fixing. 

The divergence of perturbative expansions can be expressed in an unambiguous way as the 
statement about the factorial divergence of the RG invariants pi. Anything else, for instance the 
divergence of the expansions in ([[]) and/or (||), depends on the chosen RS and has thus no direct 
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physical meaning. In my view the best strategy how to proceed in such circumstances is to follow 



the suggestion of [12], i.e. to choose at each finite order N of perturbation theory a definite 



renormalization prescription and investigate the limit of finite order approximants (S), i.e. 

lim r^(/i(A0,RS(A0). (31) 
There is no problem to choose the dependences n(N) and RS(iV) in such a way that this limit is 



finite, even for factorially divergent series with asymptotically constant sign [13]. The question is 
how to do it in a reasonably unique manner. The nonperturbative power corrections should provide 
a crucial piece of information in this respect. 
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